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ABSTRACT 


The energy and momentum equations are solved analytically to predict the heat- 
transfer coefficient of small submerged spheres in film boiling. The analysis is based 
on the postulate that the rate of entropy production is maximized. In addition, the small - 
sphere theory is matched with the film -boiling correlation for flat plates to give a heat- 
transfer correlation which applies for large, as well as small, spheres. 
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FILM BOILING FROM SUBMERGED SPHERES 
by Robert C. Hendricks and Kenneth J. Baumeister 
Lewis Research Center 

SUMMARY 

The heat-transfer coefficient for pool film boiling from submerged spheres was 
analyzed. The model developed assumed a vapor dome on top of the sphere into which 
vapor flows from the thin-film vapor region which exists on the lower portion of the 
sphere. 

In addition, the theoretical result for boiling from small spheres was matched with 
the heat-transfer coefficient from a flat plate to give a heat -transfer correlation that 
applies to both small and large spheres. 


INTRODUCTION 


Film boiling is one of the major mechanisms by which heat is transferred in cryo- 
genic systems or wherever extremely high temperature differences are involved, as in 
quenching or spray cooling. Fundamental investigations into the phenomenon of film 
boiling have considered Leidenfrost boiling of liquid drops and film boiling from flat 
plates, wires, and spheres. This report is concerned with film boiling from submerged 
spheres of relatively small diameter and some of the similarities to film boiling from 
submerged cylinders. 

Film boiling from horizontal wires (cylinders) has been analyzed theoretically and 
checked experimentally by Bromley (ref. 1). His analysis gave the following expression 
for the heat -transfer coefficient: 
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and X* is the modified latent heat of vaporization 
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(All symbols are defined in appendix A. ) To model the phenomenon, Bromley assumed 
that the generated vapor flowed laminar ly in a thin vapor annulus surrounding the heated 
cylinder. This model fit the experimental data for a range in diameters from 1/4 to 
1 inch (0. 635 to 1. 91 cm), for many common fluids. 

In 1953, Banchero, et al. , (ref. 2) furthered the understanding of film boiling. Using 
liquid oxygen as a primary test fluid, they assessed the effects of wire size and pressure 
on the film -boiling process. They pointed out that Bromley’s equation is valid over a 
narrow range of wire sizes and recommended an empirical modification: 



(4) 


Here a is nearly constant and C is a constant for each fluid tested, D is the diameter, 
and F q is the property group of Bromley's equation. 

In 1961, Berenson (ref. 3) used a critical wavelength parameter, 
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to correlate film boiling from a horizontal surface: 
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Dimensionless experimental heat-transfer coefficient, h 



Figure 1. - Heat-transfer coefficient in film boiling from horizontal cylinders as function of characteristic diameter ratio. Data from Breen and Westwater (ref. 4). Each type 
of data point represents a set of test conditions within the following ranges: Fluids, water, helium, oxygen, pentane, nitrogen, ethanol, benzene, Freon-113, isopropanol, 
or carbon tetrachloride; temperature, 72° to 2600° F (301 to 1703 K) : heat-transfer rate, 17 to 2120 Btu per hour per square foot per °F (97 to 12 200 W/(m 2 )(K)) ; cylinder 
diameter, 0.00053 to 4.6 centimeters; critical wavelength, 0.068 to 0.660 inch (0. 173 to 1.67 cm); length -diameter ratio, 0.135 to 50. 
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(b) Fluid, ethyl alcohol; wire diameter, 0.02 inch (0.0507 cm); ratio of critical wavelength to diameter, 20. 
Figure 2. - Effect of diameter on film boiling from horizontal wires. 


A year later, in 1962, Breen and Westwater (ref. 4) made another significant contri- 
bution to the understanding of film boiling. Introducing the critical wavelength x , 
they correlated all available film -boiling heat-transfer data for horizontal cylinders by 
means of a simple empirical equation. The correlation is illustrated in figure 1. This 
again demonstrates the Bromley equation to be valid over a narrow range of wire sizes. 

Baumeister and Hamill (ref. 5) developed a theoretical model for film boiling from 
small-diameter wires which followed the trend in the experimental data and gave ex- 
cellent agreement with the Breen and Westwater correlation. 

As seen in figure 1, below an l /D of 0. 1 the heat-transfer data are independent of 
wire size, whereas Bromley’s equation predicts ever decreasing values of the heat- 
transfer coefficient. The heat- transfer coefficient in this regime was correlated by 
Breen and Westwater (ref. 4) using Berenson’s results (ref. 3) for film boiling off an 
upward-facing flat plate. This correlation appears as a horizontal line in figure 1. 

Above an I /D of 2, the heat-transfer data rise at a much quicker rate than the 
Bromley equation would predict. The increase in the data slope for small wires 
(large I / D) results primarily from a change in the flow profiles around the wire, as 
displayed in the photographs in figure 2. For large -diameter wires (fig. 2(a)), the flow 
is circumferential and pseudoturbulent in nature. The flow follows a vertical pattern 
upward and leaves the wire in a large dispersed vapor column. However, for small- 
diameter wires, as shown in figure 2(b), the flow is axial into vapor domes that are peri- 
odically distributed along the wire. These domes grow as additional vapor enters, even- 
tually break away from the wire, and escape as large vapor bubbles into the bulk liquid. 
The upper photographs of figure 3 illustrate how the vapor patterns change with I /D for 
cylinders ranging from 0. 0508 to 4. 76 centimeters in diameter. Note the changes in 
vapor structure and wave patterns near the transition regime. 

Frederking, et al. , (refs. 6 and 7)performed an analysis similar to Bromley’s (ref . 1) 
for relatively large -diameter spheres. The results of this analysis did not fit their data. 
They then suggested an empirical correlation based upon the parameters of turbulent free 
convection with the property group to the 1/3 power, which fit most of their data. 

Photographs of film boiling from spheres in liquid nitrogen show a phenomenon simi- 
lar to that observed on wires of corresponding diameter (see fig. 3). Here, the small- 
diameter wires and spheres (large l /~D) both exhibit the characteristic vapor dome. As 
evidenced by the series of domes rising from the small-diameter spheres, the process 
of vapor removal is periodic in nature. Motion pictures of the phenomenon (see appen- 
dix H and film supplement) indicate that the flow around the spheres which supplies the 
growing vapor domes is laminar in nature with small capillary waves moving upward 
along the liquid -vapor interface. For large -diameter wires and spheres, however, the 
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flow appears turbulent and chaotic in nature at all positions surrounding both the wire and 
the sphere, except near the lower stagnation region. Because of this strong evidence for 
similarity between boiling on a wire and a sphere, we felt that the theoretical techniques 
used to analyze the wire could be applied to film boiling off spheres. 


Cell wavelength, A -> 
r 



(a) Small-diameter wires (UD> 2), reference 5. 



(b) Small-diameter spheres. 

Figure 4. - Film-boiling models for small-diameter wires 
and spheres. 


The theoretical model for small -diameter wires (ref. 5) was a modification of 
Berenson's original model introduced for film boiling from flat surfaces. A schematic 
of the model is shown in figure 4(a). A further modification of the vapor dome model to 
the case of film boiling from small-diameter spheres is depicted in figure 4(b). Here, 
the vapor is generated in the thin annular region at the bottom of the sphere. It then flows 
in a laminar manner into the escape dome. In the physical situation, the vapor domes 
grow as additional vapor enters, eventually breaking away from the sphere and escaping 
as large vapor bubbles into the bulk liquid. The proposed model, to be discussed in de- 
tail in the next section, assumes the existence of a time -average configuration where all 
the velocity, pressure, and temperature fields are at steady state. 


7 



Heat-transfer expressions for both large- and small-diameter spheres are developed 
herein. The data are shown to be bounded by the slip and nonslip boundary conditions; 
however, the shifts from one condition to another are unexplainable without further anal- 
ysis. For small spheres, a theoretical expression based on the model depicted in fig- 
ure 4(b) is derived and compared with the available film -boiling data. For large -diameter 
spheres, the film -boiling heat-transfer coefficient is correlated by matching the small- 
sphere and flat -plate solutions; a similar technique is used for film boiling off of a large 
cylinder. In this case, the heat -transfer coefficient is independent of geometry and de- 
pends only on a fluid property group raised to the 1/4 power. 

In addition, a motion-picture study of film boiling was undertaken to assess the na- 
ture of the vapor flow from spheres, cylinders, and a vertical flat plate. A motion- 
picture supplement C-263 has been prepared and is available on loan. A request card 
and a description of the film are included at the back of this report. 


BASIC MODEL AND EQUATIONS 

Consider the single -dome model of a submerged sphere in film boiling, as illustrated 
in figure 5. The physical situation could be a solid metal sphere supported on a wire 



Figure 5. - Schematic model of film boiling off submerged sphere, 
single-dome case. 
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and immersed in liquid nitrogen - provided, of course, that the proper temperature dif- 
ference is maintained. The material evaporated at the interface flows about the sphere 
towards the vapor reservoir. The reservoir begins at the point designated as (R Q + 5, 9*) 
in figure 5. 

The single -dome model appears to represent film boiling from small spheres (low 
Bond number, Bo) as shown in figure 3. However, for large -diameter spheres (large 
Bond numbers) the single dome appears to give way to multiple domes, such as occur off 
a horizontal plate facing upward. Apparently, as with boiling off horizontal wires, as 
the size of the sphere beco mes l arger than the critical wavelength A , from hydrodynamic 
stability theory ^A c = 2n-D/|/Bo^, the single dome tends towards multiple domes. To 
maintain a tractable solution for the large -sphere case, the conduction area under the 
single dome is assumed to be replaced by multiple -dome film boiling. This multiple - 
dome film -boiling regime, for large spheres, is assumed to be similar to film boiling 
from an upward-facing flat plate. 

There are, of course, several possible modes of heat transfer as the sphere tempera- 
ture tends toward the Leidenfrost temperature (ref. 8). These mechanisms (nucleate 
boiling, etc. ) are beyond the scope of our investigation; hence, our model is restricted 
to the developed film -boiling regime. 

The following assumptions are made in developing the single- or multiple -dome 
models: 

(1) Rotational and vibration effects due to the sphere are small. 

(2) Surface capillary waves, while influencing the boundary, do not significantly 
alter the heat transfer or the laminar nature of the flow. 

(3) The vapor gap thickness is constant. 

(4) The model has complete symmetry with respect to the coordinate. 

(5) Because of the low Reynolds number, the flow of vapor is laminar and incom- 
pressible; and the inertia and body force terms in the Navier -Stokes equations, as ap- 
plied to the vapor, are neglected. Justification for this assumption can be found in 
reference 9. 

(6) Radiation is negligible. 

(7) The velocity and temperature profiles are in steady state. 

(8) At any instant of time, the sphere is at an average temperature T , and the 
evaporating liquid is at the saturation temperature T . The properties of the flow field 
are evaluated at the film temperature 


T f = 


T + T 
w s 


and are constant (an assumption that has worked quite well, ref. 10). 
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(9) The convective terms in the energy equation are neglected; however, the latent 
heat of vaporization is corrected to accommodate this assumption. The correction is 
usually of the form 


X* = X 



( 6 ) 


Such assumptions are apparently valid, based on the work in references 8 to 11. These 
authors indicate the major mode of heat transport to be conduction; when the convection 
term is retained, the solution yields a factor analogous to the effective latent heat of 
vaporization X*. Thus, it is assumed that, in the energy equation, 


Xe9T <V r f^T « j<_ _3_ I 2 3T\ 
r 36 r 3r 3r \ dr) 


(7) 


2 . 
r sin 


i- sin e 9T 


36 


36 


2 3r 


2 3T' 
3r 


(8) 


The surrounding liquid is at the saturation temperature. Thus, all the heat reaching 
the liquid produces vapor. 

The governing equations are as follows (ref. 12): 

(1) Momentum: 
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(2) Energy: 


(3) Continuity: 



0 = i- i- (r 2 V ) + — 1- (V fl sin 6) 
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The boundary conditions are separated into two categories, which represent the 
single -vapor -dome and multiple -vapor -dome regimes, respectively. 


Regime 

Surface 

Single -dome model 

Multiple -dome model 

Equation 

number 

Thin film 

ll 

o 

v r =0, v e = 0 

T = T 

w 

(14) 


r = R q + 5 

e < e* 

V r = V r (6) ; V e = 

0, T = T s 

(15) 


e = o 

R 0 - r s R o + 6 

v e =o 

(16) 


e = e* 

r = R q + 6 

P=P d 

(H) 

Vapor dome 

s* < e < ir 

Heat transfer by conduction 

^ a (Ra' VBo d ) 1//4 

(18) 


where 0* is defined (fig. 4) as the transition from the thin film to the vapor dome, 5 
is the vapor gap thickness, and is the dome pressure. 

The boundary conditions are incomplete at this point since 6, 9*, and V (5) are 
unknowns. Hence, three additional mathematical constraints are necessary to make the 
problem tractable. 

(4) Static pressure balance: 


One additional constraint requires the sphere to be in static equilibrium. Small 
vibrations of the sphere are neglected. The static equilibrium condition requires that 
the average pressure P in the vapor at the liquid-vapor interface resulting from the 
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weight of the supported liquid, surface tension, and atmospheric pressure balance the 
average pressure due to flow (see appendix C). Summing the forces acting at the inter- 
face gives 



The total pressure in the vapor gap P can be calculated from the momentum equations, 
and P is established in appendix B, along with a solution to the momentum equations. 

(5) Interface energy balance: 

The second additional constraint is the interface energy balance. Because the sup- 
porting fluid is assumed to be at the saturation temperature (assumption 8), all the heat 
reaching the interface produces evaporation of the fluid material. Mathematically, this 
constraint is expressed as 


-pAV r (R 0 + 6 , 0) = 
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( 20 ) 


where A is the latent heat of vaporization and 



r=R o+ 5 


is the conduction heat flux of the boundary of the surrounding liquid. 

(6) Entropy production: 

The third constraint is that the rate of entropy production be a maximum. As 9 * 
is the parameter which controls the rate of exchange of energy (maximization of entropy 
production), 


oh 
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Equation (21) will lead to the determination of 6* as a function of Bond number Bo. 

For a constant area, h represents the rate of exchange of energy per degree change in 
temperature, which in a loose thermodynamic sence, " quasi -equilibrium, ” describes 
entropy production. 

This optimization technique has been used successfully in predicting the heat-transfer 
rates in film boiling from flat plates (ref. 13) and wires (refs. 5 and 14). 

(7) Slip and nonslip boundary condition: 

The analysis herein depends on the nonslip boundary condition, V Q | =0, equa- 

°K=<P 

tion (15). However, the analysis using an apparent or virtual slip boundary condition 

= 0 ( 22 ) 
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is carried out in appendix F. In subsequent material unless otherwise specified, the 
nonslip analysis is used. 


ANALYTICAL RESULTS 

The governing equations, along with the appropriate boundary conditions and con- 
straints, are solved in appendixes B, C, D, and E. Only the results are presented in 
the main text. 

Normally, the experimental data represent an average of more than one mode of heat 
transfer from the sphere. Thus, the theoretical heat -transfer coefficient is formulated 
as the sum of the two major modes of heat transport, film boiling and conduction through 
the vapor dome. With this formulation, theory and data can be compared. For a given 
fluid, the terms small Bond number and small sphere are used interchangeably, although 
the precise meaning of each is distinct. 

Small-Bond-Number Case 

For a small-diameter sphere, such as shown in the lower far right photograph of 
figure 3, the model shown in figure 5 applies. The heat -transfer coefficient may be 
written as 
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Film -boiling Conduction 

= (transport across + (transport through 
thin vapor film vapor dome) 

( e < e*)) 

where is the lower surface area of the sphere (6 < 0*), and h^ is the associated 
film -boiling heat -transfer coefficient for this area. In the second term, represents 
the surface area beneath the dome, and h ^ represents the conduction heat-transfer 
coefficient associated with this area. 

In order to compare theoretical and experimental results, h a n represents an 
area -weighted average heat -transfer coefficient from the sphere. Here, A^. represents 
the total surface area of the sphere, the sum of A^ and Ag. 

The solution for h gma jj is given by equation (B65) of appendix B as 


Nu' = 1/4 


-2Ra’G(Bo) 


1/4 -1 

+ y (1 + cos 0*)csc 9* 


where 


Nu’ = 



p(p - p)gD 3 A* 

Ra’ = L - 

kp AT 


(24) 


(25) 


(26) 


and Nu’ and Ra’ are modified Nusselt and Rayleigh numbers, respectively. The 
function G(Bo) (fig. 6(a)) depends on the Bond number, 


Bo = 


(p z - p)gR 0 


Buoyancy fo rces 
Surface-tension forces 


(27) 


For a fluid at its normal boiling point, Bond number varies, essentially, as the 
square of sphere radius (i. e. , as the size of the sphere). Other conditions, such as 
approaching the critical point or large variations in g, would effect a variation in Bond 
number for a fixed radius. 
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In addition, the angle 9 * in equations (15), (17), and (18) is a function of Bond 
number. The functional relation is shown in figure 6(b). As shown in this figure, the 
analysis predicts an upper bound in the value of 6 * to be 37t/4, as illustrated in fig- 
ure 6(b). The limit appears to agree with physical observations on boiling from small 
spheres (see fig. 3); however, it is also apparent from the other spheres of figure 3 
that the magnitude of Bo is quite a significant factor. 


Bond Number Range for Small-Sphere Model 

Considering the expression for the heat -transfer coefficient for small spheres, 
there remains one intriguing question. Over what range of Bond numbers can this analy- 
sis be applied? As the Bond number increases, experiments have demonstrated that a 
single dome no longer appears, rather multiple domes appear (lower photographs, 
fig. 3). 

From hydrodynamic stability theory, the most dangerous wavelength A m( j is de- 
fined as 




^P 7 - p) g\ 1/2 




(28) 


For the case of a heavy liquid over a lighter fluid, the most dangerous wavelength repre- 
sents the wavelength of a small sinusodial unstable disturbance, which has the fastest 
amplitude growth rate of all possible unstable disturbances. For the purpose of this 
analysis, it is assumed that if the circumference is less than the most dangerous wave- 
length, that is, if 


R 9 * < A . 
o md 


(29) 


the model and heat-transfer equation (24) would apply. Noting that 


^o = VBo 

A„ 271- 


(30) 


it follows that for 


Bo* < 
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the analysis is valid. ( 9 „„ can be found from fig. 6(b).) This Bond number Bo* pre- 

ITlcLX 

diets with reasonable accuracy the transition from single -dome film boiling to multiple - 
dome film boiling, as shown in figure 3. Consequently, to determine the heat transfer 
for Bo > Bo*, the analytic model must be modified. 


Large-Bond-Number (Large Sphere) Modification 


When the Bond number is greater than that given by equation (31), multiple domes 
replace the single domes (see figs. 7 and 3). Of course, for these larger Bond numbers, 



r Single-dome equivalent 



Figure 7. - Film boiling off submerged large sphere, 
multiple-dome case. 


the expression for the heat-transfer coefficient given by equation (24) would no longer 
apply. 

Since, at the present time, all available experimental data exist for Bond numbers 
near and outside the range of equation (24), the small-sphere analysis cannot be readily 
compared with data. It is desirable for practical use, however, to modify the theoretical 
analysis so that data can be correlated over a more complete spectrum of Bond numbers. 

The solution determined to within a constant a is found in appendix E. The solution 
for the large -Bond-number case (large spheres) becomes 

I fl2G(Bon 1//4 + Q 21 ^ + CQS g*) (32) 

(Ra-Vi ^) 1/4 4 LsV^J 
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Bond number, Bo 
(b) Optimization function GiBol. 

Figure 8. - Angular coordinate 6” and optimization function G(Bo) as function of Bond number. 








General Solution 


The solutions to the large- and small -Bond-number regimes (large and small 
spheres) in essence represent two asymptotic cases. In appendix E, the two solutions 
are matched at R Q 0*/^ C = 1, and a linear combination of the two solutions produces a 
general solution applicable to both the small- and the large -Bond-number regimes 


Nu . 1 

-2G(Bo) 

1/4 
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0 177 i csc 

(W t/Bo d ) 1/<4 4 

3VBo D 
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(1 + cos B *) 


(33) 


where Bo^ is Bond number based on sphere diameter to facilitate comparison to data. 

The functions G(Bo) and 6* are extended to the large Bo regime and plotted in 
figures 8(a) and (b). 

As shown in equation (33), if Bo^ « 1, the first term dominates, and the expres- 
sion for the Nusselt number reduces to equation (24), the expression for the heat-transfer 
coefficient for small spheres. For large Bond numbers, equation (33) approaches 


Nu' = 0. 177 (Ra'-\/Bo D 


\l/4 


(1 + cos 8*) + — 
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-2G(Bo)Ra’ 


1/4 


(34) 


The first term of this expression is the classical expression of film boiling from a flat 
plate facing upward times an area correction factor, while the second term represents 
film boiling from the lower portion of the sphere. 

The constant coefficients of equation (34) depend not only on the physical process but 
also on the boundary conditions. The slip and nonslip conditions are the most commonly 
considered variation in boundary conditions for this type of problem and are discussed 
in the next section. 


Slip Modification 


In the analysis presented in appendix F, it is assumed that the nonslip boundary con- 
dition vj = 0 (eq. (15)) occurs at the liquid-vapor interface. Bromley (ref. 1), 


e 


Z=cp 


however, showed that the slip boundary condition, along with the nonslip boundary condi- 
tion, bracketed the experimental data. For this particular problem, the slip boundary 
condition (eq. (15)) becomes 
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The details of the .modification are given in appendix F. The results of this modifi- 
cation are as follows: 

(1) Small -diameter sphere: Equation (24) becomes 
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-2Ra'G(Bo) 


1/4 -1 
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(2) Large -diameter sphere: Equation (32) becomes 
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(3) Combined solution: Equation (33) becomes 
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(35) 
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In the next section, the results for both the slip and nonslip cases (refs. 7, 15, 
and 16) are compared with available experimental data. 


DISCUSSION 

In reference 7, a transient cooldown technique was used to determine the film- 
boiling heat -transfer coefficient. The data were obtained from two spheres submerged 
in two fluids. The sphere diameters were 1/4 and 3/8 inch (0. 635 and 0. 95 cm) and 
the fluids were nitrogen and helium. The results of these tests are replotted in figure 9 
in terms of Nusselt number Nu and modified Rayleigh number Ra' . The helium data and 
the nitrogen data for the l/4-inch (0. 635-cm-) diameter sphere follow the l/4-power 
rule; however, note that the helium data are near the theoretical ’’slip" line given by 
equation (35), while the nitrogen data (1/4-in. - (0. 635-cm-) diam sphere) are near the 
theoretical line given by equation (24). Note that the nitrogen data for the 3/8-inch- 
(0. 9 5 -cm-) diameter sphere tend to follow a 1/3 -power rule, on the average. 
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Figure 9. - Transient film boiling off spheres, 
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Results from references 7 and 8. 



Figure 10. - Variation in Nusselt number with modified Rayleigh number for transient 
cooldown of 1-inch- (2. 54-cm-) diameter sphere in liquid nitrogen (ref. 8). 


More recently, transient cooldown data were obtained in reference 8 for liquid nitro- 
gen with 1/2 3/4-, and 1-inch- (1. 27-, 1. 90-, and 2. 54-cm-) diameter submerged 
spheres. These data indicate that several power rules may apply, depending on the mode 
of heat transfer. This is illustrated for the 1-inch- (2. 54-cm-) diameter sphere in 
figure 10. The circular symbols appear to follow the theoretical ’’slip" line (eq. (35)). 
However, much of the data tend to follow a 0. 515 or square-root power rule, followed 
by the usual large increase in heat transfer at nucleate boiling. In a private communi- 
cation, Lyle Gordon Rhea indicated that the shift to the 0. 515-power rule was a result 
of nucleate boiling on the stem which supported the sphere. 

With these limitations in mind (cf. large temperature differences and film boiling 
only) the data of reference 2 were replotted in figure 11 along with data from figure 10. 
These data represent temperature differences greater than 210° F (~117 K) and no 
change to the square-root rule. These data follow the 1/4-power rule; however, they 
tend to shift between the slip and nonslip theoretical lines predicted by equations (35) 
and (24), indicating a change in the model. 

Results similar to those reported in references 7 and 15 were obtained earlier in 
reference 16 (analyzed in ref. 6) for a 1-inch- (2. 54-cm-) diameter sphere submerged 
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Figure 11. - Transient film-boiling data for 1/2-, 3/4-, and 1-inch- (1. 27-, 1.91-, 
and 2. 54-cm-l diameter spheres submerged in liquid nitrogen. Selected data 
from reference 15; temperature difference between heater wall and saturation, 
>210° R and prior to - in the data of reference 15. 


in liquid nitrogen. From the results of that analysis, it would appear that the data again 
follow the 1/3 -power rule; however, these data might also follow the 1/4 -power rule if 
similar temperature limitations were applied. 

The data of figure 9 appear to be bounded by the slip and nonslip equations (eqs. (35) 
and (24)). Would these data be bounded for larger or smaller spheres? And what is the 
mechanism which causes the transition from one curve to the other? The answers to 
these questions appear to be based on geometry and fluid properties, which suggests 
that some critical dome size is involved. As indicated earlier, the Bond number equa- 
tion (eq. (27)) constitutes a parameter which incorporates these effects. 

Bond number is a natural parameter of this analysis and has been successfully ap- 
plied to film boiling from a wire (ref. 5); it is also natural to assume that film boiling 
from a sphere should show similar trends. The photographs of reference 7 and this re- 
port seem to agree with such a classification. The critical wavelengths for nitrogen, 
helium, water, and ethyl alcohol are tabulated in the following table (from ref. 17): 
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fluid 

Critical wavelength, 


in. 

cm 

Nitrogen 

0. 269 

0. 674 

Helium 

. 068 

. 173 

Water 

.660 

1. 68 

Ethyl alcohol 

. 379 

. 962 


Figure 3 and the film supplement C-263 available on loan and described at the back 
of this report and in appendix H illustrate how the vapor flow patterns change as a func- 
tion of diameter. Sphere size is one of the parameters comprising Bond number 

(P 7 - P)R? g 
Bo = —L 

ag c 

Thus, for a given fluid, an increase in sphere size increases Bond number. It is evident 
from the motion picture (film supplement), and also figure 3, that as Bond number is 
increased beyond A c /D = 1, surface capillary waves appear, signaling the onset of tran- 
sition. In the transition region, the flow is pseudolaminar (see appendix G). As the 
waves become undamped, multiple domes appear over the surface (fig. 3, i /D = 0. 19). 

It is clear that the single -dome boundary condition breaks down and the single -dome 
conduction region is replaced by a multidome film -boiling regime. This change leads to 
the large -Bond-number or multiple-dome model as sketched in figure 7. Because of the 
randomness of the bubbles in this regime and the magnitude of the Rayleigh number, the 
heat transfer may be considered turbulent; however, as shown in appendix D, Bond num- 
ber enters as a parameter, and thus the laminar equations can be used to predict heat 
transfer (see also refs. 4 and 5). 

The data of figure 5 are replotted as a function of l/Bo^ in figure 12. The theore- 
tical line of equation (33) encompasses the helium data and part of the nitrogen data. The 
nitrogen data for 1-inch- (2. 54 -cm-) diameter spheres (ref. 16) lies between 7 and 
28 percent above the theoretical line. However, major data trends apparently follow, as 
predicted by equation (33), and it would seem that Bo is a valid parameter. 

The data, other than that for helium, also tend to follow the empirical 1/3 -power 
rule, which may or may not be valid for larger- or smaller -diameter spheres than tested. 
The arguments which lead to the 1/3 -power rule were based on the lack of an effect of 
geometry on the data (ref. 6). At large Bond numbers, the effect of geometry should be 
small, as predicted by equation (33); moreover, the data of figure 11 indicate that the 
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Figure 12. - Experimental and predicted results for film boiling off submerged spheres. 


1/4-power rule is valid for high -temperature -difference film boiling. Thus, equation (33) 
should be valid over a wide range of submerged-sphere film boiling. 


CONCLUSIONS 

Laminar film -boiling heat transfer on a sphere has been analyzed subject to the optimum 
entropy production constraint for a large range of Bond numbers. The following results 
were obtained: 

1. The heat -transfer coefficients for film boiling from submerged spheres, as pre- 
dicted by equation (33), are in fair agreement with the experimental data. Equation (33) 
represents a linear combination of the low -Bond -number case (single vapor domes) and 
the high-Bond-number case (multiple vapor domes). For a specific fluid at its normal 
boiling point, Bond number varies as the square of the sphere radius. 

2. No comparison of theory and experiment could be made at very low Bond numbers 
because of the lack of data. Thus, more experimental data, particularly in the low- 
Bond-number regime (small spheres) are required before a conclusion can be made as 

to the validity of the models leading to equation (33). Experimental data at elevated pres- 
sures near the critical point and/or at high- or low -gravity conditions would also prove 
useful. 

3. The nonperiodic nature of the interface at the higher Bond numbers is clearly 
evident in the film supplement. This randomness raised the question as to whether the 
interface was turbulent or laminar. As the heat-transfer coefficient can be calculated by 
equation (33), heat transport by the vapor is apparently laminar. Thus, in a turbulent- 
looking process, the governing mode of heat transport may still be laminar in nature. 
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4. In view of the interface motion at the higher Bond numbers, the interface velocity 
boundary condition was changed from nonslip to slip. The data appeared to be bounded 
between the theoretical nonslip and slip lines; however, the shifts in level within these 
bounds were unexplainable except through the models leading to equation (33). 

5. A motion -picture supplement is offered which compares film boiling from plates, 
vertical and horizontal, cylinders and spheres. 

Before a conclusion can be made as to the validity of the analysis at very small Bond 
numbers, much more film -boiling data on spheres is required. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 11, 1969, 

129-01-05-17-22. 
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APPENDIX A 


SYMBOLS 


A 

a 

Bo 

Bo D 

C 

3 C 

C P 

C l,2,3,4 

Co 

Co , 
opt 

D 

E 2 

F 

F o 

F 1 

F 2 

G(Bo) 

G R 

g 

H* 

h 

h 

k 

L o 


surface area 
constant (see eq. (4)) 

Bond number based on radius, g(p^ - p)!* 2 /^^ 

Bond number based on diameter, 4 Bo 

constant, eq. (4) 

constant defined by eq. (C4) 

specific heat of vapor 

constants 

integration constant, eq. (B62) 
optimization constant 
diameter 

operator defined by eq. (B15) 

function defined by eq. (B32) 

function defined by eq. (4) 

function defined by eq. (B65) 

function defined by eq. (B18) 

optimization function, eq. (B64) 

function defined by eq. (B34), p(p - p) gD 3 /p 2 

acceleration of local gravity 

gravitational constant in Newton's law of motion 

/ 0. 34C AT 

modified latent heat of vaporization, X (1 + 2 

\ A , 

heat-transfer coefficient 
nondimensional component of force 
thermal conductivity 
transition length 


, see fig. 1 
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dome size parameter, 

Nusselt number 

Nusselt number based on radius of sphere 

modified Nusselt number, eq. (B56), 2(Nu R -1) 

n o 

total pressure 
pressure in vapor dome 
environmental pressure 
Prandtl number 
dimensionless pressure 
radius of vapor dome 
radius of sphere 
Rayleigh number, G^Pr 

modified Rayleigh number, Ra 

C AT 
P 

Reynolds number 

Reynolds number (appendix G), ^* u P 

M 

radial coordinate 

property group defined by eq. (B45) 
temperature 

temperature difference between heater wall and saturation 

sphere temperature 

film temperature, (T + T g )/2 

saturation temperature of liquid 

dimensionless velocity defined by eq. (B4) 

velocity at transition, laminar to turbulent 

radial velocity component 

theta velocity component 
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V 

w 

w* 

y* 

a 

/9 

y 

A 

6 

? 

0 

e 

e* 

K 

X 

X* 


dimensionless velocity 
factor defined by eq. (B36) 

nondimensionalizing factor defined by eq. (B37) 

length parameter 

matching parameter 

parameter defined by eq. (Dll) 

coefficient of dome size, see fig. 5 and eq. (B50) 

dimensionless gap thickness defined by eq. (B25) 

gap thickness 

dimensionless radial coordinate, see eq. (Bl) 

dimensionless temperature defined by eq. (B39) 
angular coordinate, see fig. 4 

angular coordinate at vapor dome, see fig. 4 

thermal diffusivity 

latent heat of vaporization 

modified latent heat of vaporization (X + 0. 5C p AT) 
modified latent heat of vaporization (X + 0. 4C p AT) 


critical wavelength, 



[X viscosity 

p kinematic viscosity, ju/p 

p density of vapor 

p ^ density of liquid 

a surface tension 

$ angular coordinate, see fig. 5 

<p dimensionless function defined by eq. (B29) 

p stream function 

Subscripts: 

cond conduction 
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D diameter 

d dome 

exp experimental 

fb film boiling 

fp flat plate 

l liquid 

max maximum 

md most dangerous 

opt optimum 

r radial direction 

t total 

6 6 direction 

1 film boiling 

2 film boiling from a flat plate 
Superscripts: 

— average value 

' derivative with respect to independent variable 
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APPENDIX B 


METHODOF SOLUTION 

In this appendix, the analytical solutions for the thickness of the thin vapor film and 
the temperatures, velocity, and pressure profiles within this film are determined. 
These parameters govern the heat transfer across the thin vapor film. A flow chart of 
the solution appears as figure 13, and frequent reference to this figure will prove fruit- 
ful. 



[Energy equations ^ 

[Nondimensionalization 


L 

o 

GO 

& 

< 


Area -weighted 

r 


i : 


film boiling 



[Area -weighted conduction! 
■ [Single-dome optimization | -— I 


[Multidome optimization!— Modified Nusselt number 
1 ■ 1 (large sphere) 


Modified Nusselt number 
(small sphere) 


Determine 

eigenvalue 


Modified 

Nusselt 

number 


Area -weighted multi - 
dome film boiling 


Figure 13. - Solution flow chart. 


Momentum Equations 


The governing equations (eqs. (9) to (13)) are made nondimensional to generalize the 
solution by selecting the following parameters: 



R 


o 


(Bl) 
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< 

</~S 

II 

C 1 < 

* I 4 

(B2) 

V 0 

v 0 = — 
u* 

(B3) 

2,1 o 

1 os 

ii 

* 

3 

(B4) 

pu*R 

R - ° = 1 

(B5) 



P- P 

(B6) 


*2 

pu* 


Substituting these parameters into the momentum equations (eqs. (9) and (10)) gives 


where 


0 = - — + V 2 v*, - — v r - — — ^ - — v cot 0 

^ ^ ^2 de ^2 y 


(B7) 


0 = -lgP + v 2 

? 00 ^2 de ^2 sin 2 g 


(B8) 


V 2 = _ L _ L L2 1 j_ / sin0 A. 

C 2 9 ? \ 9i 7 C 2 sin 0 00 \ de > 


(B9) 


Introducing the stream function for spherical coordinates (ref. 12, p. 131) gives 
for the velocity distribution 


v 


r 


1 

2 

r sin 0 


dxjj 

00 


(BIO) 
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V = 
0 


di{/ 


r sin 6 3r 


(Bll) 


or, in terms of the dimensionless velocities, 


v r = - 

5 ,2 


My) 


£ sin 9 


3 9 


(B12) 


V 9 


dMS,0) 


£ sin 9 3? 


(B13) 


Substituting equations (B12) and (B13) into the momentum equations (eqs. (B8) 
and (B7)) and combining to eliminate the pressure terms gives 


E*0//) = 0 


(B14) 


where 


E 2 _ d__ + sin 9 5 


1 3 


3?" 


c 


2 30 \sin 0 30; 


(B15) 


The similarity transform 


MS, 9) = F 2 (?)(l - cos 0) 


(B 16 ) 


is used to convert equation (B14) into an ordinary differential equation. This type of 
transformation has been used for solving for flow around a sphere (ref. 17, p. 217). 
Substituting equation (B16) into equation (B14) yields 


Fg(?) = 0 


(B17) 


This is a linear homogeneous equation of fourth order, the solution of which is 


f 2 (?) = c 1 + c 2 ? + c 3 ? 2 + c 4 ? 3 


(B18) 
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Therefore, 



v 


e 


T 






(B19) 


(B20) 


The constants in the previous equations must be determined from the boundary con- 
ditions (eqs. (14) and (15)). In terms of the dimensionless variables, these boundary 
conditions become 


v 0 (l, 6) = o 

(B21) 

v 0 (l + A,0) = 0 

(B22) 

< 

V 

CD 

II 

O 

(B23) 

Vj»( 1 + A,0) = V A 

(B24) 


where 


A 



(B25) 


Applying these conditions to equations (B19) and (B20) allows the four constants to be 
expressed in terms of one constant of integration. 


C 1 = -±(1 - 3 c p) 
1 2 

(B26) 

C 2 ~ ^^4 *P 

(B27) 

3 C, 


c 3 = - — - (1 + <p) 

(B28) 
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where 


cp = 1 + A (B29) 

The dimensionless gap thickness A is assumed at this time to be an unknown con- 
stant. However, as is shown later in this appendix, this requires that the interface 
energy balance be satisfied on an average over the heating surface rather than at every 
point. 

The constant can be found by satisfying the static force balance constraint. 

First, however, the pressure distribution in the vapor gap must be found. Substituting 
equations (B19) and (B20) into equations (B7) and (B8) and solving for the pressure dis- 
tribution gives 


2C 

p(£, 0) = — 3 - 6C 4 In £ + 12C 4 In (sec + 3 C (B30) 

where C„ is related to C. by equation (B28). 

0 O 

The constants C 4 and °C are determined from the static pressure balance in 
appendix C in terms of the parameters 0* and Bo; thus, for cp = 1, 



GrF( 0*, Bo) 
3X2 6 


(B31) 


Ml sin - (1 - cos 0* + 2 sin 0*) 

F(0*, Bo) = - B ° sitl — — (cos 0* 

In cos 2 — + sin 3 — 


1 ) 


(B32) 


3 C=p 


sin 0*1 


8 \Bo sin 0* p 


l 


- 2C 3 + 6C 4 In cp - 12C 4 In 



(B33) 


where 


P(P Z - p)gD 3 _ 6w 


2 7TW* 

P 


(B34) 
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(B35) 


g(p 7 - p)R„ 

Bo = i 

CTg c 

W = — (p ; - p)R^ (B36) 

3 g c 

pu* 2 R 2 

w* = 9. (B37) 

g c 

Therefore, substituting the constant in C^, Cg, and Cg, the velocity and pres- 
sure distributions are known relations in cp and 0*, with Bo as a parameter, where 
cp = 1 + A. 

Next, the energy equation is considered. Afterwards, the solutions to the energy and 
momentum equations are combined in the interface energy balance to obtain a solution for 
the heat -transfer coefficient in terms of 6*, with Bo as a parameter. 

Energy Equation 

The nondimensional form of the energy equation (eq. 12) becomes 

— k 2 — \ = 0 (B38) 

a? \ ac/ 

where 

T - T 

0 = — (B39) 


Equation (B39) can be solved directly and, when the boundary conditions (eqs. (14) 
and (15)) are used, yields 


0 = — (l - (B40) 

The temperature gradient at each surface becomes 
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= £ 
A 


(B41) 


30 

3 ? 


l£=l 


30 1 

3? £=1+A 


(B42) 


Interface Energy Balance 

The velocity and temperature distributions have been expressed, up to this point, 
in terms of an unknown parameter, the dimensionless vapor gap thickness A. The 
interface energy balance (eq. (20)) is now used to determine the value of this parameter. 

First, the radial velocity at the interface must be determined. Substituting the 
values of the constants to Cg into equation (B19) and evaluating the velocity at the 
interface (£ = 1 + A) gives 

v r (cp, 0) _ v A _ i A _ 3( p [ 3 
C 4 " C 4 " 2 V cp 2 

where is a function of 6* and the parameter Bo (see eq. (B31)). 

Nondimensionalizing equation (20) results in the interface energy balance 


(B43) 


v 


A 


30 

3? 


?=1+A 


(B44) 


where 


S 


i r c P ( T w- T s) 

Pr A* 


(B45) 


and A has been replaced by A* (see assumption 9 and eq. (3)). 

Substituting the expression for v^ from equation (B43) into equation (B44) along 
with the expression for the gradient (eq. (B42) and solving for A with cp « 1, gives 


A 


4 


2cpS 

C 4 (0*,Bo) 



3 

2RaF(0*, Bo) 


(B46) 
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The problem has been reduced to one of determining 0*. 

Now the analytical solution for the total heat transfer through the thin vapor film and 
the single vapor dome is developed. The heat transport through the vapor dome and the 
thin vapor film (film boiling) will be area -weighted, and the maximization principle will 
provide 9* as a function of Bond number for the low -Bond -number case. 

The heat -transfer coefficient and Nusselt number can be constructed from the mo- 
mentum and energy solutions with 6* as a parameter. 

The heat -transfer coefficient to the sphere is defined as 

(B47) 


h(T w - T s ) = -k £i| 


or, in nondimensional form, 



hR 

o 


30 


k 



Substituting equation (B41) into equation (B48) gives 


Nu- 


R, 


=£= 1+1 


(B48) 


(B49) 


This equation represents the Nusselt number for heat conduction in the thin vapor layer 
beneath the sphere (0 < 0*). 


Approximation for Heat Conduction Through Upper Sphere 

The vapor dome may be roughly approximated by concentric spheres of radius Ro 
and (R q + y R q <p sin 0*), as shown in figure 5. The approximation for A, shown by the 
dashed line in figure 5, is 


A = R 0 y<p sin 9* 


(B50) 


where y is a constant to be discussed in detail in appendix E. Substituting equation (B50) 
into the right side of equation (B49) gives 
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( N x) 


R, 


= 1 + 


cond ^ Sin 


(B51) 


The average Nusselt number for the sphere follows as the linear sum of the area -weighted 
film -boiling and conduction Nusselt numbers. 


0 


Nu 


R 


NUr 


v o A. 


Jfb 


( N X) 


Me > 9*) 

cond 


(B52) 


Jcond 


The fraction of the surface area in film boiling for a given 9* can be found directly by 
integration: 


x = 


r / 

o M 


■9* 


dA 


1 - cos 9* 


/ 2n ri\ 

•S(\ 


(B53) 


dA 


' 0 •'O 

The conduction area then is, of course, 

A co n d< e>8 *) = (' - x > A t 


(B54) 


When equations (B49), (B51), and (B53) are substituted into equation (B52), the average 
Nusselt number becomes 


Nu 


R. 


/i - cos e*\ cp 

_L 

A + cos 0*\ 4 , 1 \ 

Li 2 AJ 

fb 

\ 2 ) \ ytp sin 9*) 


(B55) 


Jcond 


Rearranging yields 


Nu’ 


= 2 ("Nu l) 1 - cos 9* + (1 + cos g*)csc 6 s1 
' R o A cpy 


(B56) 


Substituting the value of A from equation (B46) and assuming cp « 1 yields the following 
form of equation (B56): 
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Nu' = 1 _._ ? s e * F(e*,Bo)J 1 / 4 + (1 + cos 0*)y -1 CSC 0* 


(B57) 


This expression is a function of 9* and Bo. As Bo is dependent on the type of evapora- 
tive fluid and geometry, in any given experiment, the unknowns are reduced to 6* alone. 

To determine 9*, we apply the optimization criterion to Nu' . For the case of a 
fixed geometry and constant properties, this criterion may be written as 


D e * Nu ' = 0 


D;* Nu’ < 0 


(B58) 


Differentiating Nu' and equating to zero gives the differential equation 

i-l 


4 D e* ln 


-F(0*, Bo) 


sin 6* 

1 - cos 6* 


1 r_p Ra t "I 1/4 

- ZKa F(0*,Bo) 

4 L 3 J 


(1 - cos 0*)y 


= 0 


(B59) 


Note that the last term in equation (B59) may be rewritten by using equation (B46) as 


(1 - cos 8 *) y 


(B60) 


However, we have assumed A to be small, so equation (B59) reduces to 


D e* ln 


-F(«*,Bo) 


4 sin 9 * _ q 

1 - cos 0* 


(B61) 


where 9* is not equal to or greater than n . Rather than complete the above differentia- 
tion, equation (B61) can be integrated directly to give 

F(1 - cos 0*) 4 = Co (B62) 

which gives a value for 9* that maximizes the Nusselt number according to equa- 
tion (B58). Here Co is a constant of integration which represents the locus of boundary 
conditions determined by the interrelation of 9* and Bo. 

Substituting equation (B62) back into equation (B57) gives 
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(B63) 


Nil' = - 
4 



+ (1 + cos 0*)y 


-1 


esc 9 * 


The optimum (or maximum) Nusselt number now depends on Co, 9 *, and Ra'. The 
maximum values of Co can be determined by s/ solving equation (B62) and plotting Co as 
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Figure 14. - Families of solutions to F(1 - cos O 4 ')^ * Co for various bond 
numbers and angular coordinate 9 :5 values. 


a function of 0* for various Bond numbers Bo at a given Ra' . This is illustrated in 
figure 14, with the dashed line denoting the locus of Co maximums for various Bond 
numbers at some Ra’. The optimum integration constant Co and 9 * values for 
various Bond numbers are given in figure 6. The locus of Co ^ is labeled G(Bo) and 
is required in the solution of Nu’. 


Co opt - G(Bo) 


(B64) 


Based on these results, the average Nusselt number for the single -dome case becomes 


Nu’ = - 
4 


-2Ra’G(Bo) 


1 1/ 4 


+ F^Bo) 


(B65) 


where G(Bo) and F^(Bo) = (1 + cos 9 *) y - * esc 9 * are determined from figure 6 for a 
given Bo. The effect of y is discussed in appendix E. 
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APPENDIX C 


STATIC PRESSURE BALANCE 

The pressure distribution within the vapor gap is obtained directly from the solu- 
tion to the momentum equations and is repeated here for convenience. 

p(C, 0) = — 3 - 6C 4 In e + 12C 4 In /sec |J + 3 C (B30) 


where 


3C 

Co — {1 + <p) (B28) 

J 2 

and 3 C and C 4 are, as of this point, undetermined constants. In this appendix these 
constants are determined, using the boundary condition (eq. (17)) and the force balance 
(eq. (19)), with 0* and <p as parameters. Figure 15 illustrates the contributing pres- 
sure terms. 
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The pressure in the vapor dome at the exit surface 9* is, in dimensional terms, 


p - p + §£_ | PS p 

d ° R g d 
K d g c 


where 


= R Q cp sin 6* 


(Cl) 


(C2) 


The dome pressure represents the sum of the static pressure P , the surface- 
tension pressure difference (2a/R Q ), and the vapor head. Therefore, the constant 
in equation (B30) may be evaluated as follows: 


p (cp,9*) = — — = — - - 6C, In cp + 12C. In (sec — \ + 
pu. 2 v * \ 2/ 


(C3) 


or 


C = P o + 


2a 


Pg R o<P \ 2C„ 

+ sin 0* — 2 ± 


R 0 <P sin 9* g c 


pu* 


2 cp 


+ 6C^ In cp - 12C^ In |sec — ^ 


(C4) 


And substituting equation (C4) into equation (B30) yields 


p (? , 9) = 2 C 3 (l-l)- 


6C 4 In 


<P. 


+ 12C 4 In 


sec — 


sec 


9* 


+ 


P o + 


2a 


pu* 


R cp sin Q* g 
o^ °c 


+ PI Rtt, sine*) (C5) 


As one of the constraints, the average pressures across the interface must be bal- 
anced. Therefore, to evaluate equation (19), we need to establish P, the average pres- 
sure exerted by the liquid upon the interface. The pressure at any point (£, 0) becomes 
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P = P + P 1 R<5lS + (R + 5)(p. - p)&- (cos 6 - cos 9*) + 2cr (C6) 

e c «C R 0 + 5 

Total pressure P represents the static head P , the liquid head to the base of the 
vapor dome p R^g/g^ the difference in liquid and vapor head to any point below (cp, 9*), 
and the pressure drop due to surface tension, which must be balanced by the pressure 
within the vapor gap. Thus, equation (C6) can be integrated over the sphere to obtain 

P = /p dA 
'/dA 


re* 

l < R o 


+ 6) 2 sin 9 


P +p ,R n — + (R,, + 5 )(p 7 - p) — (cos 9 - cos 0 *) h — — — — 
° L g 1 g R +5, 

& C & C O -I 

'0* 


f 6 2 

/ (R + 5) sin 9 d0 


For a constant gap thickness and a fixed fluid, equation (C7) becomes 


d0 

- (C7) 
(C8) 


p = p g + R d P z 


S- + - 2ct - + 

S c R o + 5 


(R o + 6)(p^ -p)j- 

fc> 


(1 - cos 9*) 


(C9) 


Integrating the pressure (eq. (C5)) over the sphere gives 
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But by equation (19) 


fiz* 

J g„ 


<p sin 0 de d<f> = Ij 


(C12) 


Substituting equations (C9) and (Cll) into equation (C12) and integrating gives 


12C. [ In cos 2 — + sin 2 — 

4 \ 2 2 j 


cos 0* - 1 

*2 
p u* 


2a 


R Q <p sin 0* 


(1 - sin 9*) 


S r 


- (pj - P) R n ^ sin 0* - - (p - p) Rep (1 - cos 0*) (C13) 

1 Sr, 2 1 g ° 


Introducing the Bond number Bo 


Bo - 


(p z - p)gR^ 


(B35) 


the weight of the displaced fluid (sphere alone) 


w = H(p, -p)R 3 £- 
3 1 6c 


(B36) 


w* = 


* 2„2 
pu* R q 

S„ 


(B7) 


into equation (C13) and solving for yields 


= -J ?.?.. (cos 9* - 1) 
% 32ttw* 


4(1 - sin 9*) 
2 

Bo cp sin 9* 


(1 - cos 0* + 2 sin 0*) 


i 2 0* - 2 0* 

In cos — + sin — 


(C14) 
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(Cl 5) 


C = _W£_ F ( Bo< ^ Q *) _ 

4 32jrw* 


G R F(0*,Bo<p 2 ) 

3X2 6 


Thus, remains a function of Bo, 0*, and the gap thickness cp = 1 + A to be determined 
by other conditions. 
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APPENDIX D 


LARGE-SPHERE MODIFICATION 


In this appendix, the large -Bond-number (large sphere) case is considered. The 
single dome is replaced by a multiple -dome region, and the basic assumptions are de- 
lineated. Then, the overall heat transfer through the thin film and the multiple -dome 
region is determined. The two regions are area -weighted, and the maximization princi- 
ple is applied to give 0* as a function of Bond number and a matching parameter. An 
equation for heat transfer from large -Bond -number spheres is presented in the form 


h 

n large ~ 



+ 



Film boiling 
= (transport across 
thin vapor film, 
part 1) 


Film boiling from upward- 
+ facing flat plate (transport 
through multiple -dome 
region) 


(Dl) 


where a is an eigenvalue or matching parameter. 

Film boiling from a large sphere is modeled in figure 7. In solving the large -sphere 
regime for the optimum heat-transfer coefficient (or, equivalently, for the Nusselt num- 
ber), the following assumptions are made: 

(1) Flow from the lower stagnation region to the multiple -dome region is described 
as in the single -dome regime. The solutions of the momentum and energy equations are 
assumed to be valid to the point of optimization of the heat -transfer coefficient which 
determines 0* = 0*(Bo). 

(2) The conduction region of the single -dome model is replaced by a multiple -dome 
film -boiling regime. The heat-transfer coefficient in this region is proportional to that 
for a flat plate facing upward; the determination of the constant of proportionality in 
optimizing the heat-transfer coefficient leads to an eigenvalue problem. 

(3) The heat transfer may be expressed as a area-weighted sum of the heat transfer 
in the lower film -boiling region (0 < 0* in fig. 7) and film boiling off the upper portion 
of the sphere (0 > 0* as described in assumptions (1) and (4)): 


Nu’ 

2 


Nu x 


-±+ NUr 


*t 



(D2) 
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where the Nusselt number for boiling off the lower portion of the sphere is 


Nu' = Nu - 2 


with 


NUj 



and 



1 - cos 9* 
2 


(D3) 


(D4) 


(D5) 


For the top of the sphere, 


and 


Nu 0 = a 


t i ~\ 1/4 

■ BA' VBo ] 


77 


/ 


(D6) 



1 + cos 9* 
2 


(D7) 


and a represents an eigenvalue or matching parameter, to be determined in matching 
the small- and large-sphere regions, as (R0*/X )— 1. Note, a factor of 1/2 was intro- 
duced in the left side of equation (D2) so that the Nusselt number would be based on 
diameter (see eqs. (B55) and (B56)). 

(4) The ratio of vapor gap thickness to sphere radius A is assumed to be sufficiently 
small that cp = 1 + A ~ 1. At this point, it is helpful to consult figure 13 to determine 
the cause of subsequent events. Equation (D2) may be rewritten as 


Nu’ = 1 T C . 9 S -A— + (1 + cos 9 *) 
A 




a 


(D8) 
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The dimensionless gap thickness A may be written as (see eq. (B46)) 


1_ 

A 


_4 

2Sy 


1/4 


wF(fl*, Bo) 
2xSx32ttw* 


, 11/4 _ i 
4 


2F(0*,Bo) 


L 3 


Vbo J 


1/4 


(ita’ VBo) 


1/4 


(D9) 


\l/4 


Substituting equation (D9) into equation (D8) and factoring out ^Ra’ Vbo ) gives 

1 1/4 /o\ 1/4 


Nu' 


1 - cos 0* 


(Ra' Vbo) 


1/4 


2F(0*,Bo) 


L 3 


Vbo J 


£ (1 + cos 0*) 
4 


(DIO) 


where the modified matching parameter (3 is defined as 

/> - 4 (^) V4 - V 

V2/ 1/4 


(Dll) 


Equation (DIO) will now be optimized with respect to e* to give the optimum large- 
sphere heat-transfer coefficient (equivalently Nusselt number). Differentiating equa- 
tion (DIO) with respect to 9* yields the characteristic differential equation: 


D 0 * ln|~F(0*,Bo)| 


4 sin 6* 

1 - cos 0* 


r 

_ 

l/4"j 

•si - A 

Vbo 

1 

V 

-F(0*,Bo) 

J 


In appendix E, the matching parameter is shown to be 


a 


IT 


1/4 


0. 21 


(D12) 


(E4) 


Applying equations (E4) and (Dll) to equation (DIO) gives an equation applicable to the 
large -Bond -number case (large sphere) 



1 

4 


- 2G(Bo) 
3 Vbo 


1 1/4 


+ 0. 21(1 + cos 0*) 


(D13) 
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APPENDIX E 

COMBINED SOLUTION 


In appendix B, an analytical solution for film boiling from spheres was found for the 
low-Bond-number case. In essence it represents an asymptotic case applicable only 
where the single vapor dome forms over part of the sphere. In appendix D, a semi- 
analytic solution for film boiling from spheres was determined to within a constant for 
the large -Bond-number case. 

Small -bond -number case: 


Nu’ = - 
4 


— Ra’G(Bo) 
3 


1/4 -1 

+ (1 + cos 9*)y CSC 0* 


Large -bond-number case: 


/ Nu' \ 1//4 _ 1 

-2G(Bo) 

\Ra' VBo/ 4 

_ 3 Vbo _ 



(1 + cos 0*) 


(B65) 


(El) 


where a is a parameter to be determined. 

To determine a Nusselt number applicable to both the large- and small-sphere re- 
gimes (large and small Bo), a must be determined such that the solutions match in the 
transition regime (film supplement and appendix G). Now the matching parameter a is 
determined by requiring the large- and small-Bond-number solutions to match at R Q 9*/ 
A. = 1. Superimposing the small-sphere solution for 0* as a function of Bo (fig. 7) 
onto figure 16 indicates that the two solutions are compatible; that is, 


R 0* 
o 


when the modified matching parameter j3 becomes 


and from equation (Dll) 


13= 0. 93 


- a .= /£\/2\ 1//4 ~ o.21 

1/4 \4/\3/ 

71 \ ' 


(E2) 


(E3) 


(E4) 
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A single equation for film boiling off spheres is obtained by forming a linear combi- 
nation of the large- and small -Bond -number cases (large and small spheres). Linearly 
combining equations (B65) and (El) and substituting equation (E4) for gives an 

equation applicable to both large and small Bond numbers 


/ Nu- y/ 4 _ i 

-2G(Bo) 

1/4 

0 21 + csc ^ 

\Ra' VBo / 4 

3 Vbo 

i 

a/ 4 

(Ra’ VboJ 


(1 + cos 9*) 


(E5) 


where G(Bo) and 6* = 0*(Bo) are determined as functions of Bo from figure 8. 

In order to compare the results of this analysis to that of the cylinder and experi- 
mental data, it is most expedient to convert equation (E5), in part, to Bond number 
based on diameter rather than on radius. 


Bo = 



(E6) 


Substituting equation (E6) into equation (E5), except for the function G(Bo), results in 


^Eta' VBOq) 


Nu ' _ _ 1 
1/4 4 


-2G(Bo) 


Vbo 


D. 


1/4 


CSC 0* 


y 


-l 


0. 177 + 


, 1/8 


Ra' 


VBo 


D 


1/4 


(i + cos e*) 


(E7) 


The expressions for the Nusselt numbers require the evaluation of the constant y 
for both the small and large spheres (eq. (B65)) and implicitly for the large sphere (see 
fig. 7 for single -dome equivalent and fig. 5 for A). The parameter y, in effect, de- 
termines the average thickness of the vapor dome through which heat is transferred by 
conduction. This parameter (y /4 ' ) which appears in equation (E7) was assumed to 
be 1, in agreement with the very large single vapor domes shown in figure 3. With this 
assumption, equations (B65), (E5), and (E7) become equations (24), (32), and (33), re- 
spectively. 
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APPENDIX F 


SLIP CASE 

When fluid motion at the interface is present, the problem becomes quite complex; 
however, when the problem is limited to the single -dome model with slip at the inter- 
face, a tractable model results. The boundary conditions are modified as follows: 

dvJcp,9) 

— H .= 0 

3? 


v 0 (l,e) =0 L 
Vj*(i, e) = o 


(-F1) 


v^{cp,e) = v A 


As the solution of the governing equations is the same, only the values of the constants, 
Cj, Cg, Cg, and change. Substituting equation (FI) into equation (B18) and solving 
for the constants in terms of yields 


C 4 2 

Ci = — i(l - 3 cp 2 ) 

1 2 

C 2 = 3cp 2 C 4 

C 3 = -!c 4 (1 + cp 2 ) 

2 

Using these constants to solve for the interface velocity gives 

C 4 A 3 

v. = (3 cp + 1) 

A p 

2 cp 


(F2) 

(F3) 

(F4) 


(F5) 


This, of course, alters the magnitude of the interface energy balance over that pre- 
viously determined for the nonslip case. Using equations (F5) and (B44), the interface 
energy balance for the slip case is 
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V 


A 



S 

cp A 



(3 <p + 1) 


Solving equations (F6) and (F7) for A 

A 4 = (s\ 

For the situation where cp is approximately 1, 



(F6) 


(F7) 


(F8) 


(F9) 


where is identical with the previous definition of (eq. (B31)). A comparison of 
equations (F9) and (B46) lead to a Nusselt number ratio for the slip and nonslip cases: 


Nu 


slip 


Nu 


nonslip 


nonslip' 
'slip 


1/4 


= V2 


(F10) 


Therefore, the slip condition (eq. (B65)) may be written as 


Nu’ = 




- - G(Bo)Ra’ 

- 3 


1/4 


+ F i (Bo) 


(Fll) 
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APPENDIX G 


SURFACE WAVES 

Hsu (ref. 18) studied film boiling off a vertical flat plate. In order to predict local 
heat -transfer data, a transition length L q had to be determined. Hsu assumed that the 
transition from laminar to turbulent flow occurred at a characteristic Reynolds number 
of Re* = 100. Knowing Re*, the transition length L q may be calculated as 

L = R e * 

° 2kAT 

where 


y* = 


2/i 2 Re* 
gp(P z - P) 


1/4 


The variation of L Q with temperature difference is illustrated in figure 17 for two fluids, 
water and nitrogen. While both fluids may possess minimums, that of nitrogen is quite 
pronounced near 200 K. This also suggests that the surface becomes more stable near 



Figure 17. - Critical height for transition to turbulent flow. 
Theory from reference 3; pressure, 1 atmosphere. 
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I I 


Wave amplitude growth 

and bubble domes in Liquid nitrogen 

turbulent region % 





- Capillary waves in 
laminar region 


(a) On surface of vertical flat plate. 




(c) On surface of sphere. 

Figure 18. - Illustration of flow profiles. 


the Leidenfrost temperature, an observation verified by experiment. Flow on a vertical 
plate is illustrated in figure 18(a). Note the capillary waves in the laminar region which 
roll off the bottom of the plate. These waves are most easily demonstrated by the film 
supplement (available on request). 

Some more recent observations of thin films (ref. 19) indicate the existence of 

several regimes - laminar, pseudolaminar, transition, pseudoturbulent, and turbulent. 

From these observations, it is evident that Re* = 300 could represent an upper bound 

4/3 

to the transition regime or an increase in L Q by 3 ' . 

Thus, for film boiling of nitrogen off a vertical flat plate at a temperature difference 
of 480 K, L q would fall in the range 1. 02 < L q < 4. 4 cm (see fig. 17). Within this 
length and for some distance beyond, the surface capillary waves and disturbances will 
not strongly affect the heat transfer (i. e. , it remains laminar in nature). 

Analogous results are found in film boiling off cylinders, as seen in figure 18(b) and 
the film supplement. In this case, the heat transfer is predicted in reference 5 with a 
laminar analysis and by optimizing the heat transfer with respect to a critical wavelength. 
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Transition 
123° 



Laminar and 
pseudolaminar 


(a) Sphere diameter, 3/8 inch (0.95 cm). 



(b) Sphere diameter, 1/2 inch (1. 27 cm). 

Figure 19. - Laminar, pseudolaminar, and 
transition regimes for spheres submerged 
in nitrogen. 


Even though the surface appears to be quite turbulent (fig. 2), the heat transport is gov- 
erned by the laminar equations and the Bond number. 

The sphere is no exception. Capillary waves travel over the surface (see fig. 18(c) 
and the film supplement). At large Bond numbers, the surface appears to be quite turbu- 
lent (fig. 3). Pursuing the concept that R Q 0 and L Q are directly related, we can pos- 
tulate the laminar, pseudolaminar, and transition regimes for spheres submerged in 
nitrogen (fig. 19). 

A comparison of figure 19 and figure 3 indicates some capillary wave activity in the 
pseudolaminar region, but the waves do not appear to grow. In the transition regimes, 
however, these waves appear to be growing, and vapor is sometimes sheared off, which 
forms bubbles. 

It is to be concluded that while a process may appear turbulent, the governing charac- 
ter of the heat transport may still be laminar. 



APPENDIX H 


MOTION- PICTURE SUPPLEMENT - A DESCRIPTION 

The motion-picture supplement ^resents visual observations on film boiling of liquid 
nitrogen from three different geometries: the flat plate, oriented in the horizontal and 
vertical positions, a series of horizontal cylinders of different diameters, and a series of 
spheres of various diameters. The diameters of the cylinders and the spheres are 
similar to permit comparison of the boiling phenomenon from large and small diameters. 
The interaction of diameter and interface stability, Bond number, is discussed and il- 
lustrated. By changing the size of the cylinder and the sphere, the viewer can visualize 
the effects of Bond number, see also figure 3. For a more complete discussion of the 
stability effects, see appendix G. 


Effect of Surface Disturbances 

For each of these geometries, the relation between small disturbances and a critical 
wavelength is discussed to determine if the disturbances grow or decay (fig. 20). A sta- 
bility map which summarizes some essential results of hydrodynamic stability theory is 
presented as figure 21. If the disturbance wavelength is less than the critical wavelength, 
the disturbance will decay; however, if the disturbance wavelength is greater than the 
critical wavelength, it will grow, and an unstable interface results. The wavelength 
for maximum -growth is predicted and labeled the most dangerous wavelength ^ mc j; how- 
ever, experimental data tend to fall to the right of A as seen on figure 21. 


Liquid 


-Wavelength, A- 

I 9 

Surface- 

tension 


Surface 

disturbance 



Vapor 


Figure 20. - Surface disturbances related to critical 


A < A, 


disturbance is damped; 
, disturbance grows; where A c = ^ 


wavelength „ c 
A > A, 


<P i - pig 
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R 2 Rj 



Figure 22. - Gaussian curvature 
of typical surface element. 


Vertical Flat Plate 

A surface with no curvature (fig. 22), where — °° and Rg — is placed in liquid 
nitrogen, and the ensuing flow phenomena are photographed. A thin vapor layer quickly 
forms over the surface and begins to grow. At the first transition point, some instabi- 
lities are noted, while at the second transition point, Hsu (ref. 18) notes a significant 
change in heat-transfer coefficient. This signifies a change to turbulent heat transfer. 

The observation is quite important because even though the interface is unstable, the 
heat transport can still be governed by the laminar equations; however, the significant 
shift in the coefficient at the second transition is irrevocable evidence of a change in 
mechanism. 

Unfortunately, the plate was not instrumented and only the observations of film 
boiling are available. 


Horizontal Flat Plate 

While this surface has no curvature, its orientation with respect to the gravity vec- 
tor changes, which affects the interface geometry. In this case, the interface may be con- 
sidered to be composed of many vapor domes, as depicted in figure 23. The heat vapor- 
izes the liquid at the interface, and it flows into the vapor domes. These vapor domes, in 
reality, grow and are released into the fluid reservoir. This action is seen in the motion 
picture. The multiple domes reach an optimum spacing according to the critical wave- 
length for the fluid (ref. 13). A periodicity of the vapor domes can be seen in the film. 


rVapor dome 
Liquid 1 9 . 



Heating surface 


Figure 23. - Model of film boiling from 
horizontal surface. 
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Horizontal Cylinder 


Here is finite, however R 2 — 00 . This represents a significant change in geo- 
metry. A secondary geometric change is invoked as R^ — A c /2; in this regime the 
basic flow about the cylinder is altered. 

Small cylinder. - When R.. « A /2, the vapor flows parallel to the axis of the 
cylinder and escapes in a regular periodic manner through vapor domes, as seen in 
figure 24(a). 


K r vapordomes 

.r-y 


Vapor flow pattern 



(a) Small. 



-Vapor flow 
pattern 


Figure 24. - Vapor flow patterns for large and 
small heating cylinders. 


Large cylinder . - When R^ > A c /2, the vapor flows around the cylinder and escapes 
in a chaotic manner near the upper portion of the cylinder (see fig. 24(b)). 

As shown in the main text, the terms large and small are relative to a given fluid; 
it is the Bond number that determines the regime. 

To illustrate the changes in vapor patterns, several different size cylinders were 
placed in liquid nitrogen and photographed. 


Sphere 

For this geometry, R^ = R^. This represents another significant change in geometry. 
It is apparent from the motion picture and the models presented in this report that there 
exist large and small sphere regimes analogous to that found for the cylinder. Again, 
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large and small are relative terms, and it is the Bond number which determines the 
regime. 

Small spheres . - This regime is characterized by a single vapor dome, « A c /2, 
forming and departing in a periodic manner. The phenomenon is modeled in figure 5. 

The motion picture indicates some surface capillary waves even for the single -dome 
cases, which, for the larger spheres, grow and lead to the multiple -dome regime. 

Large spheres . - Many vapor domes appear on the upper portion of the sphere, and 
the periodicity is dubious (R^ > X c / 2); however, the model (fig. 7) assumes some regu- 
larity to form a tractable problem. This was also done for the large cylinder. 

To illustrate how the flow patterns change, several different-size spheres were 
placed in liquid nitrogen and photographed. No heat-transfer data were taken. 

In the closing sequence of the film, the flow patterns encompassing the "large" 
sphere and cylinder are compared; at the same time, the patterns about the "small" 
sphere and cylinders are compared. Cross comparisons and comparisons between the 
geometries can also be made. 
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